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ABSTRACT 



Aims. We aim to understand the formation of dense cores by magnetosonic waves in regions where the thermal to magnetic pressure 
ratio is small. Because of the low-ionisation fraction in molecular clouds, neutral and charged particles are weakly coupled. Ambipolar 
diffusion then plays an important role in the formation process. 

Methods. A quiescent, uniform plasma is perturbed by a fast-mode wave. Using 2D numerical simulations, we follow the evolution 
of the fast-mode wave. The simulations are done with a multifiuid, adaptive mesh refinement MHD code. 

Results. Initial perturbations with wavelengths that are 2 orders of magnitude larger than the dissipation length are strongly affected 
by the ion-neutral drift. Only in situations where there are large variations in the magnetic field corresponding to a highly turbulent 
gas can fast-mode waves generate dense cores. This means that, in most cores, no substructure can be produced. However, Core D of 
TMC-1 is an exception to this case. Due to its atypically high ionisation fraction, waves with wavelengths up to 3 orders of magnitude 
greater than the dissipation length can be present. Such waves are only weakly affected by ambipolar diffusion and can produce dense 
substructure without large wave-amplitudes. Our results also explain the observed transition from Alfvenic turbulent motion at large 
scales to subsonic motions at the level of dense cores. 

Key words. MHD - Shock waves - ISM: clouds - ISM: TMC-1 - Stars: formation 
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1. Introduction 

Observations show that molecular clouds are highly structured 
(e.g. Blitz & Stark 1986). Furthermore, emission-line profiles of 
molecular tracers such as CO and CS are considerably broader 
than their thermal line widths indicating the presence of highly 
turbulent motions (e.g. Falgarone & Phillips 1990 ). Since molec- 
ular clouds are threaded by magnetic fields, it is natural to sup- 
pose that these observed line widths are due to magnetohydro- 
dynamic (MHD) waves (Arons & Max |1975i >. 

Large-scale, three-dimensional (3D) simulations of turbulent 
gas motions (e.g. Ballesteros-Paredes & Mac Low 2002; Padoan 
& Nordlund |20Tj21 Gamrnie et al. 120031 Li et al. 120041 Galvan- 
Madrid et al. 120071 ) show that dense cores with statistical prop- 
erties similar to those observed can indeed be formed in this 
way. By following the evolution of a single MHD wave, Falle 
& Hartquist (120021 in ID and Van Loo et al. (120061 ) in 2D found 
that dense cores could be generated by the excitation of slow- 
mode waves. This process works on different length-scales and 
can also explain the formation of substructure within cores them- 
self (Van Loo et al. 120071 . 

In all these simulations, the plasma is treated as a single 
ideal plasma in which the plasma is perfectly coupled with the 
magnetic field. However, the low ionisation fraction in molecu- 
lar clouds (Elmegreen 119791 ) implies that the plasma and mag- 
netic field are actually weakly coupled on the scale of the cores. 
The charged particles then drift through the neutral particles giv- 
ing rise to ambipolar diffusion (Mestel & Spitzer [T9561 ). Large- 
scale, 3D simulations including ambipolar resistivity (Padoan et 
al. 120001 Oishi & Mac Low 120061 ) show that dense cores still 
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arise in those conditions. Lim et al. ( 120051 ) investigated in ID the 
effect of ambipolar resistivity on the evolution of a single MHD 
wave. They found that ambipolar diffusion affects the evolution 
of waves with wavelengths up to about a thousand times the dis- 
sipation length-scale and thus has a significant effect on much of 
the observed structure in star formating regions. 

In this paper we extend the model of Lim et al. ( 120051 ) to two 
dimensions. The governing equations and the initial condition 
are given in Sect. [2] We then examine the density structure in 
the numerical calculations (Sect. [3). and discuss the relevance of 
these results in Sect. [4] 

2. The model 

2.1. Governing equations 

Since the ionisation fraction within molecular clouds is low, the 
plasma needs to be treated as a multicomponent fluid consisting 
of neutrals, and charged particles. Here we assume the charged 
particles to be ions and electrons only. In the limit of small mass 
densities for the charged fluids, their inertia can be neglected. 
The governing equations for the neutral fluid are given by 



dp„y„ 
dt 



% + V(p„v„) = 0, 
ot 

+ V(p„V„V„ + p n ) = J X B, 



(1) 



with V = (d/dx,d/dy). For the charged fluid j, the equations 
reduce to 



dt 



— + Vipjvj) = s jn , 



(2) 



ajPj(E + vj x B) +pjp„Kj„(\„ - \j) = 0. 



2 
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where aj is the charge to mass ratio, Kj„ the collision coefficient 
of the charged fluid j with the neutrals, sj„ the mass transfer rate 
between the charged fluid j and the neutral fluid, B the magnetic 
field, E the electric field and J the current given by J = £j ffjPjVj. 
We adopt an isothermal equation of state p = pT for all fluids 
with T the isothermal temperature. Note that the units are chosen 
so that all factors of An and c disappear from the equations. 

Within translucent clumps the ionisation fraction, X,-, 
changes as a function of the neutral density (e.g. Ruffle et 
al. |1998l >. For a visual extinction of Ay = 3, the ionisation frac- 
tion decreases from X/ = 1CT 4 for «// < 1000 cm -3 to = 10~ 7 
for «// = 5 x 10 4 cm -3 , roughly following a 1 /«// law. As the ion- 
isation and recombination rates are sufficiently high in molecular 
clouds, ionisation equilibrium is quickly reached. We therefore 
adopt mass transfer rates Sj„ so that the fractional abundances of 
ions and electrons vary as l/«//. 

The Hall parameters of the ions and electrons, i.e. the ra- 
tio between their gyrofrequency and the collision frequency 
with the neutrals, are large in most astrophysical situations 
(Wardle[[998). The ions and electrons thus gyrate many times 
around a magnetic field line before they collide with a neutral 
particle. These collisions affect the evolution of the magnetic 
field which is governed by 



((V x B) x B) x B 



(3) 



<9B / 
— -Vx(v„xB) = Vxk 

where r a is the ambipolar resistivity (Falle 2003). As the elec- 
trons have a Hall parameter which is larger than the ions, the 
ambipolar resistivity is given by 

B 2 



PiPnKin 

Equations Q] - [3] are solved with an adaptive mesh refine- 
ment code using the scheme described by Falle ( 120031 1 which 
we extended to two dimensions. This scheme uses a second- 
order Godunov solver for the neutral fluid equations (Eqs. QJ 
and b). The charged fluid densities are calculated using an ex- 
plicit upwind approximation to the mass conservation equation 
(Eq.[2k), while the velocities can be calculated from the reduced 
momentum equations (Eq. Wp). The magnetic field is advanced 
explicitly, even though this implies a restriction on the stable 
time step at high numerical resolution due to the ambipolar re- 
sistivity term, i.e. At < Ax 2 /4r a (Falle |2Q03T >. 

The code uses a hierarchy of grids such that the grid spacing 
of level n is Ax/2", where Ax the grid spacing of the coarsest 
level. The solution is computed on all grids and the difference 
between the solutions on neighbouring levels is used to control 
refinement. More specifically, if the mapped down coarse cell 
value of the neutral flow variables differs by more than 10% from 
the fine cell value, the grid is refined. The loss in efficiency due 
to the time step restriction is then partially balanced by the gain 
due to adaptive mesh refinement. We also implemented the di- 
vergence cleaning algorithm of Dedner et al. (120021 1 to eliminate 
the errors due to non-zero V ■ B. 



we study the formation of dense inho- 



2.2. Initial conditions 

Like Lim et al. (120051 
mogeneities by following the evolution of a fast-mode wave in a 
uniform plasma in which the magnetic pressure dominates over 
the gas pressure. Slow-mode waves are not excited in the initial 
state. They have a very low propagation speed, and will there- 
fore remain localised near the boundary until they are generated 
inside the cloud by nonlinear steepening of the fast-mode wave. 



The quiescent, uniform background plasma conditions are 
given (in dimensionless units) by 



p„ = l,p e = 5 X 10 



10- 3 ,v„ 



0,B X = I, By 



0.25, 



and an isothermal temperature T = 0.015 for all fluids. 
Furthermore, we adopt a e = -2.0 x 10 12 and cr, = 10 8 for 
the electron and ion charge to mass ratio, and K e „ = 4.0 x 10 5 
and Kj„ = 2.0 x 10 4 for the electron and ion collision coeffi- 
cients. These values correspond to the properties of a translucent 
clump with optical extinction Ay — 3, a neutral number density 
of 500 cm -3 and a magnetic field of 30 /iG for a core tempera- 
ture of 10 K. Note that we adopt m„ = 2m// and m, = 10m// for 
the neutral and ion mass. The adopted value for the ion mass lies 
well within the observed range, i.e. from * 3m// when H^ is th e 
dominant ion to w 30m// when it is HCO + (Caselli et al. 120021 1. 

A non-linear fast magnetosonic wave propagating in the 
positive x-direction is then superposed onto the uniform back- 
ground. The initial state needs to be calculated using the method 
described in Lim et al. (120051 ) and Van Loo et al. (120071 ). since 
a simple linear approximation is not valid. This method requires 
that a wave satisfies 



dP 

dx 



oc r 



where P is the vector of the primitive variables p, v and B, 
and Tf the right fast-mode eigenvector (see Falle & Hartquist 
120021 ). By specifying the profile for v y (x, (/), all other primitive 
variables can be easily calculated. We assume that, for all flu- 
ids, the (/-component of the velocity associated with the wave 
changes sinusoidally with x. Furthermore, we introduce a sinu- 
soidal phase-shift with respect to the (/-direction to make the flow 
two-dimensional, so that 



v y (x,y)-Ai sin 



2k / 2mi 

— x + At sin 

A \ " I A 2 



where A is the wavelength of the fast-mode wave, A 2 the wave- 
length of the perturbation in the (/-direction and A 2 the amplitude 
of the phase shift in the (/-direction. A i is the amplitude of the 
(/-component of the velocity and is chosen so that the maximum 

amplitude of the total velocity in the wave ( ^jv 2 + u 2 ) is a spec- 
ified value A. We perform calculations for A = 0.7, 1.0 and 1.5. 
As the Alfven speed, v a = B/ yp, is roughly 1.0 in our simula- 
tions, these models correspond to sub-, trans- and super- Alfvenic 
perturbations. As in to Van Loo et al. (12006b . we set A 2 = 0.01/1 
and A 2 = A. The spatial variations in the fast-mode waves with 
different wavelength A are then identical and any differences in 
the dynamical evolution is due to the physical processes and not 
the initial conditions. 

Due to the ambipolar resistivity in a weakly ionised plasma, 
a fast-mode magnetosonic wave of long wavelength A dissipates 
on a time-scale 

A 2 

td = — . 

r a 

Since we are interested in the generation of dense structures 
associated with slow-mode waves, it is necessary that the ini- 
tial fast-mode wave can generate these slow-mode waves before 
it dissipates. The time for the wave to steepen into a shock is 
roughly 



A 
2A X 



(4) 
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which gives the following condition 

r a 0.025 
i>> 2A7 = — 

for dissipation to be unimportant initially. Like Lim et al. (12005 \ , 
we study the evolution of fast-mode waves with A = n and 107r 
which satisfy the above condition. The numerical domain is < 
x < A and < y < A with periodic boundary conditions. We use 
4 grid levels with a finest grid resolution 800 x 800. 

We also do simulations for A — IQOn, but these are done 
with the ideal MHD code described in Falle dl991t . This is be- 
cause we cannot resolve shocks if A = 100 n with the above 
resolution. The physical shock thickness, which is of the order 
of the ambipolar dissipation length Id, is smaller than the finest 
grid spacing. (The dissipation length U is the wavelength of the 
wave that dissipates in one wave period, i.e. Cftj — A with Cf the 
fast-magnetosonic speed and Id = 0.05 in our simulations.) This 
itself is not a problem as the numerical scheme smears out the 
shock front over a few grid cells. However, the magnetic field 
then appears discontinuous on the grid and this cannot be han- 
dled by our numerical procedure (see Falle 120031 ). Although we 
could use higher numerical resolution this increases the compu- 
tational cost considerably and is unnecessary since the results of 
Lim et al. ( 120051 1 show that dissipation is negligible for A > 30n 
and the flow can be treated as a single ideal fluid. 

3. Results 

3.1. Large wavelengths, e.g. A = IQOn 

It is useful to first consider the evolution of fast-mode waves in 
the absence of ambipolar dissipation. This is not only an appro- 
priate description for waves with wavelengths larger than ss 307r 
(see above), it also allows us to properly describe the effect of 
ambipolar dissipation. 

Van Loo et al. (2006) performed similar simulations in which 
they followed small-amplitude fast-mode waves. These simu- 
lations show that the first dense structures are formed by the 
non-linear steepening of the wave and that the subsequent in- 
teraction of the fast-mode shock with these dense regions gener- 
ates additional structure. Since non-linear steeping occurs earlier 
for larger initial amplitudes (see Eq. @), we would expect dense 
structures to appear earlier for such waves. Furthermore, it was 
shown that larger wave amplitudes result in larger density con- 
trasts, i.e. p max oc Aj. Figure \T\ shows that this dependence still 
holds for finite-amplitude waves. 

The density structure for finite-amplitude waves is more fila- 
mentary than produced by the small-amplitude waves. This is 
a result of the more rapid steepening of the fast-mode wave 
and the concurrent faster growth of the perturbation in the y- 
direction. Figure [2] shows that some substructure with p > 10, 
or Yin > 10 4 cm -3 (p = 1 corresponds to n# = 10 3 crrT 3 , see 
Sect. 12. 2\ , is present within these filaments. As these density 
structures have densities higher than the threshold for line emis- 
sion of NH3, they can, in a traditional sense, be referred to as 
dense cores (e.g. di Francesco et al. 120071 1. In our simulations, 
dense cores arise where a shock interacts with another shock or 
with an already present dense region. In such interactions slow- 
mode waves are excited which are associated with density en- 
hancements (e.g. Falle & Hartquist l20021 Van Loo et al.[2006). 

Although the initial fast-mode wave is not affected by am- 
bipolar dissipation, it can be expected that the dense cores are. 
However, slow-mode waves associated with dense cores do not 
generate large perturbations in the magnetic field (see e.g. Falle 
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Fig. 1. Temporal evolution of the maximum normalised number 
density for A = 1007T with different wave amplitudes: A = 1.5 
(solid), 1.0 (dashed) and 0.7 (dotted). The normalisation factor 
is nil - 10 3 crrT 3 . 
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Fig. 2. Density structure for A = IQQn with wave amplitude 
Ai = 1.5 after 4 wave periods. The contour lines indicate nor- 
malised number densities higher than 5, 10, 15 and 20. The same 
normalisation factor as in Fig.Q~]is used. 



& Hartquist |2002| l. (Note that this means that the magnetic field 
within a dense core is roughly the same as the background mag- 
netic field.) As ambipolar diffusion only acts on the magnetic 
field, slow-mode waves propagate without any significant damp- 
ing (Balsara 119961 ). The lifetimes of gravitationally unbound 
dense cores are, therefore, not determined by the ambipolar dis- 
sipation, but by the dispersion time-scale of a core. This disper- 
sion time scale is given by fdisp ~ f/v a , where r is the core radius, 
and is generally short compared to the wave period of the initial 
wave. Figure \T\ shows that high-density contrasts are present for 
several wave periods. Slow-mode excitation and the associated 
dense core formation only ceases when the fast-mode shock has 
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time (in wave periods) 

Fig. 3. Same as Fig.[T]but for A = lOn. 



dissipated away most of its energy. The formation of structure is 
therefore driven by the fast-mode wave. 

3.2. Intermediate wavelengths, e.g. A = lOn 

For the wave-amplitudes studied here, a fast-mode wave with 
a wavelength A — 10n steepens into a shock before ambipolar 
diffusion becomes important. Although this suggests a similar 
evolution for waves with A = 10n as for A = lOOn, there are 
some important differences. Figure[3]shows that the high-density 
structures develop later than for A = lOOn (see Fig.[T). Also, the 
maximum densities within cores are roughly a factor of 5-10 
smaller. Note that for A = lOn, A/lj - 200 n and that ambipolar 
resistivity is thus important even for waves with wavelengths that 
are rather large compared to the dissipation length. 

It is not the fast-mode wave that is significantly affected by 
ambipolar diffusion, but the initial perturbation perpendicular to 
the propagation direction is. So far, we have neglected its role 
in the dynamics. The initial ^-dependent contribution of the per- 
turbation due to the phase-shift of the fast-mode wave has the 
same wavelength as the fast-mode wave, but the amplitude of 
the variations in the magnetic field due to this phase shift are 
about two orders of magnitude (or more) smaller. The variations 
in the (/-direction will therefore dissipate more quickly than the 
those in the x-direction. Evidently, the variations perpendicular 
to the magnetic field grow more slowly. Figures |4] a and b show 
the difference in the growth of these variations for A = \00n and 
A = Wn. 

As a result, it becomes more difficult to produce large density 
perturbations. The largest density inhomogeneities form within 
dense filaments due to the interaction of these regions with 
shocks (e.g. Van Loo et al. 120061 ). The core density thus depends 
on the density of the parent filament and of the shock strength. 
In this case, not only do the filaments have smaller densities than 
for A = 1007T, but the shock strength is continuously decreasing 
due to the ambipolar dissipation. This also explains why there is 
only a short period of time (i.e. ss 3-4 wave periods) for which 
we find peak densities with p > 10 or > 10 4 ctrT 3 . 

3.3. Short wavelengths, e.g. A = n 

In the previous section, we found that ambipolar dissipation has 
a moderate effect on the evolution of waves with a wavelength 
A = lOn. For shorter wavelengths, i.e. A — n, the impact is 
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Fig. 4. The normalised density structure at half a wave period 
for A = 1007T (upper panel) and A - \Qn (lower panel) with 
A=1.5. The density range is set between 1 and 4 with the same 
normalisation factor as in Fig.Q] 



greater even though the wave is two orders of magnitude longer 
than the dissipation length. The initial ^-dependent contribu- 
tion to the perturbation is, as can be expected from the previ- 
ous case, not large enough to develop any significant variation 
in that direction. This means that there are only variations in the 
x-direction (see Fig.|5}. The evolution of the wave is then nearly 
one-dimensional. 

Ambipolar dissipation causes rapid decay of the fast-mode 
wave: the shock wave loses its energy within a few wave peri- 
ods. The generation of density inhomogeneities then stops, as 
can be seen in Fig. [6] The dissipation process, however, also 
contributes to the formation of density structures at short wave- 
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Fig. 5. The normalised density structure after 3 wave periods for 
A-7T with A = 1.5. 
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Fig. 6. Same as Fig.[TJbut for A — n. 



lengths. Lim et al. (2005) showed that, for A * n, ambipolar 
dissipation plays a major role in the excitation of slow-mode 
waves. Figure [6] indeed shows that the density contrasts arise 
earlier than for A = 10tt. However, ambipolar diffusion is still 
more efficient in reducing the peak density which does not be- 
come much higher than p = 5 or n# = 5 x 10 3 cm -3 . 

To break the one dimensionality of the flow and form cores, 
the initial phase shift needs to be larger than in our standard 
model: the variations of the magnetic field in the (/-direction 
should be of the same order as the .^-direction. Additional simu- 
lations for A — n with A 2 > 0.1 A produce dense cores with sizes 
of the order of the dissipation length. The maximum density is 
roughly 2-3 times higher than in the standard model (see Fig. [6j 
even for the largest wave amplitudes. This means that fast-mode 
waves with wavelengths that are 2 orders of magnitude longer 
than the dissipation length can only form cores in a highly tur- 
bulent plasma. 



4. Discussion and conclusions 

In this paper we studied the formation of density structures in a 
weakly ionised plasma by following the evolution of a fast-mode 
wave. This adequately describes the formation of dense cores in 
a clump which is perturbed at its edge. Slow-mode waves prop- 
agate at low speeds and remain near the boundary. Fast-mode 
waves propagate much faster and produce slow-mode waves 
within a clump either due to non-linear steepening or interac- 
tions with denser regions. However, in a weakly ionised plasma, 
the ion-neutral drift is an important dissipation process. We find 
that ambipolar resistivity strongly affects the evolution of fast- 
mode waves with wavelengths up to 2 orders of magnitude larger 
than the dissipation length. 

For an MHD wave in which the velocities are of the order of 
the Alfven speed, the ambipolar dissipation length-scale is given 
by (e.g. OsterbrockQMT] Kulsrud & Pearce [TPr59i l 



l d = 100 



B 



10" 



30fiG j\ Xj )\ tin 



1000 



3/2 



AU, 



(5) 



so that A = n » 100 Id corresponds to 0.05 pc. Note that 
this wavelength is roughly consistent with the size of a dense 
core. Fast-mode waves with A = 100/j do not produce signifi- 
cant density enhancements, even for large velocity perturbations. 
However, in a sufficiently inhomogeneous medium, these waves 
can still produce moderately-dense cores with sizes of the order 
of the dissipation length scale. These cores have a short life- 
time which is of the order of decades. The existence of such 
short-lived, small-scale clumps has been inferred in the diffuse 
interstellar medium from the time variability of interstellar ab- 
sorption lines (e.g. Dieter et al. 1 19761 Crawford l20021 >. 

Small-scale features have also been detected in the cores 
of cold dense clouds such as in Core D of TMC-1. Peng 
et al. d!998t identified, within TMC-l's Core D, 45 micro- 
clumps with sizes of 0.03-0.06 pc and masses of 0.01-0.15 Mq. 
Most of these microclumps have masses below the Jeans mass. 
Therefore, they cannot have formed by gravitational instability. 
Instead these microclumps may be generated by magnetosonic 
waves. 

With a core density of ;?h ~ 6 x 10 4 crrT 3 and fractional ion- 
isation Xj » 10' 6 (Caselli et al. [19981 and B * 15//G (Turner & 
Heiles 120061) . Id = 10 AU in Core D. The dissipation length is 
thus considerably smaller than its size of 0.1 pc, in fact waves 
with wavelengths up to 3 orders of magnitude larger than the 
dissipation length scale can be present in the core. In our mod- 
els these wavelengths correspond to A = I0n. Density pertur- 
bations are generated relatively easily by such waves, even for 
sub-Alfvenic velocity perturbations. This fits well with emission 
line observations of e.g. CS, NH3 and CO which show a substan- 
tial non-thermal broadening that is not highly supersonic with 
respect to the sound speed in H2 (e.g. Fuller & Myers [19921 ). 
Magnetosonic waves can thus generate the substructure within 
dense cores. 

However, TMC-l's Core D has an atypical fractional ioni- 
sation. Most cores in the sample of Caselli et al. (Q998 ) have 
an ionisation fraction which is 1 to 2 orders of magnitude lower 
than in Core D of TMC-1. The dissipation length then increases 
by the same order of magnitude (see Eq. [5]) making it unlikely 
that substructure is generated within those cores. 

Lim et al. (2005) showed that ambipolar resistivity only has 
a small, or even a negligible, influence on the evolution of fast- 
mode waves with A > 1000Zrf(= I0n). The slow-mode waves 
generated by the fast-mode wave are also not affected by am- 
bipolar diffusion (Balsara |1996l ). The flow can then be described 
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by a single ideal MHD fluid. A fast-mode wave with a wave- 
length similar to a clump size (1007T ~ 5 pc) produces dense 
cores with «h = 10 4 - 10 5 ctrT 3 and radii of the order 0.1 pc 
which are similar to observed properties (e.g. Jijina et al. |1999l >. 

It is clear that, with these physical properties, most dense 
cores in our simulations have masses above the thermal Jeans 
mass and are, thus, gravitationally bound. However, their masses 
are still lower than or of the same order as the magnetic Jeans 
mass. Magnetic pressure support then prevents the gravitational 
collapse of the core. Including self-gravity in our models will 
affect the dynamical evolution of the gas, but we expect it to 
be only a local effect, i.e. the density of a dense core will in- 
crease. As a consequence, the dense core will become super- 
critical and can undergo collapse. Simulations of the evolution 
of weakly ionised, magnetised, self-gravitating clouds justify 
this picture. Basu & Ciolek ( 12004| > and Li & Nakamura ( 120041 1 
in 2D and Kudoh et al. (120071 1 in 3D showed that sheet-like 
clouds fragment within a dynamical timescale to form supercrit- 
ical dense cores surrounded by subcritical envelopes. However, 
further simulations are required to confirm the effect of self- 
gravity. 

In this paper we examined the global changes to the den- 
sity structure due to ambipolar diffusion. However, our results 
also give a qualitative idea of its effect on the velocity structure. 
For typical clump properties the fast-mode and Alfven wave 
with wavelengths corresponding to dense cores dissipate their 
energy quickly (in » 10 4 yr). Slow-mode waves, however, are 
not subject to dissipation (Balsara fl 9961 ). The velocity dispersion 
on small scales is therefore dominated by the slow-mode wave 
while all types of waves contribute to the velocity dispersion at 
large scales. By examining the eigenvectors for the slow, fast 
and Alfven waves in a low-/? plasma (see e.g. Falle & Hartquist 
120021 1. it can be shown that the velocity perturbation of a fast- 
mode wave with a moderate density perturbation is of the order 
of the Alfven speed, whereas it is only of the order of the ther- 
mal sound speed for slow-mode waves. This qualitative picture 
agrees well with the observed transition from supersonic mo- 
tions in molecular clouds to subsonic motions in dense cores 
(Mvers [T983l . 
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